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Abstract 

While the naive factorization assumption works well for many two-body nonleptonic B 
meson decay modes, the recent measurement of B — > _D(*)°M° with M = tt, p and u> shows 
large deviation from this assumption. We analyze the B — > D^M decays in the pertur- 
bative QCD approach based on factorization theorem, in which both factorizable and 
nonfactorizable contributions can be calculated in the same framework. Our predictions 
for the Bauer-Stech-Wirbel parameters, |a2/ai| = 0.43±0.04 and Arg(a2/ 'a\) ~ —42° and 
1 02/01 1 = 0.47 ± 0.05 and Arg{a2/a\) ~ —41°, are consistent with the observed B — > Dir 
and B — > D*tt branching ratios, respectively. It is found that the large magnitude 1 0,2 1 
and the large relative phase between 02 and a\ come from color-suppressed nonfactor- 
izable amplitudes. Our predictions for the B° — > D*°p°, D*°lu branching ratios can be 
confronted with future experimental data. 
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1 INTRODUCTION 



Understanding nonleptonic B meson decays is crucial for testing the standard model, and 
also for uncovering the trace of new physics. The simplest case is two-body nonleptonic B 
meson decays, for which Bauer, Stech and Wirbel (BSW) proposed the naive factorization 
assumption (FA) in their pioneering work pQ. Considerable progress, including generalized 
FA (21 El H] and QCD-improved FA (QCDF) [5], has been made since this proposal. On the 
other hand, technique to analyze hard exclusive hadronic scattering was developed by Brodsky 
and Lepage [6 i based on collinear factorization theorem in perturbative QCD (PQCD). A 
modified framework based on kx factorization theorem was then given in |H] , and extended 
to exclusive B meson decays in [§J El EU EI • The infrared finiteness and gauge invariance of 
kr factorization theorem was shown explicitly in [13J. Using this so-called PQCD approach, we 
have investigated dynamics of nonleptonic B meson decays [Til Ho"! ITS] . Our observations are 
summarized as follows: 

1. FA holds approximately for charmless B meson decays, as our computation shows that 
nonfactorizable contributions are always negligible due to the cancellation between a pair 
of nonfactorizable diagrams. 

2. Penguin amplitudes are enhanced, as the PQCD formalism includes dynamics from the 



region, where the energy scale /i runs to J Am^ < m&, A = — being the B meson 



and b quark mass difference. 

3. Annihilation diagrams contribute to large short- distance strong phases through the (S + 
P) (S — P) penguin operators. 

4. The sign and the magnitude of CP asymmetries in two-body nonleptonic B meson decays 
can be calculated, and we have predicted relatively large CP asymmetries in the B — * 
K^tt HUTU and tttt modes PS [El EE]. 

All analyses involving strong dynamics suffer large theoretical uncertainties than we would 
like. How reliable are these predictions? This can be answered only by comparing more of 
our predictions with experimental data. For this purpose, we study the B — > D^M decays 
in PQCD, where M is a pseudoscalar or a vector meson. A meson is massive, and the 
energy release involved in two-body charmed decays is not so large. If predictions for these 
decays agree reasonably well with experimental data, PQCD should be more convincing for two- 
body charmless decays. Since penguin diagrams do not contribute, there are less theoretical 
ambiguities, such as the argument on chiral enhancement or dynamical enhancement. Checking 
the validity of PQCD analyses in charmed decays is then more direct. 

Note that FA is expected to break down for charmed nonleptonic B meson decays ^H]- FA 
holds for charmless decays because of the color transparency argument: contributions from the 
dominant soft region cancel between the two nonfactorizable diagrams, where the exchanged 
gluons attach the quark and the antiquark of the light meson emitted from the weak vertex. 
For charmed decays with the light meson replaced by a meson, the two nonfactorizable 
amplitudes do not cancel due to the mass difference between the two constituent quarks of 
the D" meson. Hence, nonfactorizable contributions ought to be important. This observation 
further leads to the speculation that strong phases in the B — > D^*'M decays, if there are 
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any, arise from nonfactorizable amplitudes. In charmless decays, strong phases come from 
annihilation amplitudes through the (S + P)(S — P) penguin operators, since nonfactorizable 
ones are negligible as explained above. Annihilation amplitudes should not be the source of 
strong phases for charmed decays, which do not involve the (S + P)(S — P) penguin operators. 

In this paper we shall apply the PQCD formalism to the two-body charmed decays B — > 
D^M with M = 7r, p and to. PQCD has predicted the strong phases from annihilation ampli- 
tudes for charmless decays, which are consistent with the recently measured CP asymmetries in 
the B% — ► 7T + 7r _ modes. It is then interesting to examine whether PQCD also gives the correct 
magnitude and strong phases from nonfactorizable amplitudes implied by the isospin relation 
of the B — > D(*'M decays. Compared to the work in jllj . the contributions from the twist-3 
light meson distribution amplitudes and the threshold resummation effect have been taken into 
account, and more modes analyzed. The power counting rules for charmed B meson decays, 
constructed in [20J, are employed to obtain the leading factorization formulas. It will be shown 
that nonfactorizable contributions to charmed decays are calculable in PQCD, and play an im- 
portant role in explaining the isospin relation indicated by experimental data. The predictions 
for the B — > D*°p°, D*°uj branching ratios can be confronted with future measurement. 

In Sec. II we review the progresses on the study of two-body charmed nonleptonic B meson 
decays in the literature. The PQCD analysis of the above decays is presented in Sec. Ill by 
taking the B — > Dtc modes as an example. Numerical results for all the B — > D^M branching 
ratios, and for the extracted BSW parameters a\ and a<i are collected in Sec. IV. In Sec. V 
we compare the PQCD approach to exclusive B meson decays with others in the literature. 
Sec. VI is the conclusion. Appendix A contains the explicit expressions of the factorization 
formulas. 



2 REVIEW OF PREVIOUS WORKS 

2.1 PQCD Approach to5^ L>W Form Factors 

To develop the PQCD formalism for charmed B meson decays, we have investigated the B — > 
transition form factors in the large recoil region of the meson j20j. We briefly review 
this formalism, which serves as the basis of the B — > D^M analysis. The B — > transition 
is more complicated than the B — ► 7r one, because it involves three scales: the B meson 
mass tub, the meson mass m D (*), and the heavy meson and heavy quark mass difference, 
A = niB — mf, ~ m D {*) — m c of order of the QCD scale Aqcd, t^dm (wic) being the meson 
(c quark) mass. We have postulated the hierachy of the three scales, 

m B > m D M > A , (1) 

which allows a consistent power expansion in m D (*)/mB and in A/m D (*). 

Write the B (Z>W) meson momentum Pi (P 2 ) in the light-cone coordinates as 

The picture associated with the B — > transition is shown in Fig. where the initial state 
is approximated by the bd component. The b quark decays into a c quark and a virtual W 
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Figure 1: B 
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Figure 2: Lowest-order diagrams contributing to the B — > DW form factors. Quite a bit of 
momentum must be transfered to the spectator d quark through the hard gluon exchange. 

boson, which carries the momentum q. Since the constituents are roughly on the mass shell, 
we have the invariant masses kf ~ 0(A 2 ), i — 1 and 2, where k\ (k 2 ) is the momentum of the 
spectator d quark in the B (£)W) meson. The above kinematic constraints lead to the order of 
magnitude of k\ and k 2 [20] . 



K ~ (A, A, A), 

„ (jHJL A ^ KA ) . (3) 

\m D (*) m B J 

The lowest-order diagrams contributing to the B — > form factors contain a hard gluon 
exchange between the b or c quark and the d quark as shown in Fig. |21 The d quark undergoes 
scattering in order to catch up with the c quark, forming a meson. With the parton 
momenta in Eq. (jSJ), the exchanged gluon is off-shell by 

(h - k 2 f ~ — — ^— A 2 , (4) 

which has been identified as the characteristic scale of the hard kernels. Under Eq. we have 

(*) 

mB/iTiD ^ 1) anc i the hard kernels are calculable in perturbation theory. It has been found 
that the applicability of PQCD to the B — > Z?W transition at large recoil is marginal for the 
physical masses m B and m D (») j2Ti] . 

Infrared divergences arise from higher-order corrections to Fig- 12 The soft (collinear) type of 
divergences is absorbed into the B (D^) meson wave function (p B (x 1 ,bi) (0 D (*)(x 2 , b 2 )), which 
is not calculable but universal. The impact parameter b\ (b 2 ) is conjugate to the transverse 
momentum k\x (k 2 x) carried by the d quark in the B (D") meson. It has been shown, from 
equations of motion for the relevant nonlocal matrix elements, that </\b(^i>&i) (</>£>(*) (x 2 , 62)) 
has a peak at the momentum fraction x\ = k^ / Pf ~ A/mg (x 2 = k 2 / P 2 ~ K/m D (*)) |2TT| . 
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The form factors are then expressed as the convolution of the hard kernels H with the B 
and meson wave functions in k? factorization theorem, 

F BD( *' \q 2 ) = J dx 1 dx 2 d 2 b 1 d 2 b24>B(xi,bi)H(xi,x 2 ,b 1 ,b2)(j) D (*)(x2,b2) ■ (5) 

The D^*> meson wave function contains a Sudakov factor arising from kx resummation, which 
sums the large double logarithms a s ln 2 (m£&2) to all orders. The B meson wave function also 
contains such a Sudakov factor, whose effect is negligible because a B meson is dominated by 
soft dynamics. The hard kernels involve a Sudakov factor from threshold resummation, which 
sums the large double logarithm a s In 2 x\ or a s In 2 X2 to all orders. This factor modifies the 
end-point behavior of the B and meson wave functions effectively, rendering them diminish 
faster in the small X\^ region. 



2.2 End-point Singularity and Sudakov Factor 

It has been pointed out that if evaluating Fig. 0m collinear factorization theorem, an end-point 
singularity appears |21j . In this theorem we have the lowest-order hard kernel, 

H<®{ Xl ,x 2 )<x-^, (6) 

from the left diagram in Fig. EJ which leads to a logarithmic divergence, as the meson dis- 
tribution amplitude behaves like </>£>(*) (^2) oc X2 in the small X2 region. This singularity implies 
the breakdown of collinear factorization, and kx factorization becomes more appropriate. Once 
the parton transverse momenta kj- are taken into account, Eq. (jUJ) is modified into 

4 

" (,)(l " l! ' k,T ' kir)0< M >m | + (k 1T k 2T )^ 2 m! + k iT ] ■ (7) 

A dynamical effect, the so-called Sudakov suppression, favors the configuration in which kx is 
not small [S|. The end-point singularity then disappears as explained below. 

When an electron undergoes harder scattering, it intends to radiate more photons. Hence, 
the scattering amplitude for radiating no photons must be suppressed by a factor, whose effect 
increases with the electron energy. In QED it is the well-known Sudakov suppression factor, 
an amplitude for an electron not to emit a photon in hard scattering. In the current QCD 
case of the B — > D^*' transition, it is the c-d quark-antiquark color dipole that undergoes 
hard scattering. When the color dipole is larger, in intends to radiate more gluons. Since the 
final state contains only a single ZJW meson, the real gluon emission is forbidden in the hard 
decay process. Similarly, the transition amplitude must involve a Sudakov factor, whose effect 
increases with the size of the color dipole, i.e., with the separation b between c and d. That 
is, the configuration with a smaller separation b or with a larger relative transverse momentum 
kj> is preferred in the B — > ZJW transition at large recoil. Then the virtual particles involved 
in the hard kernel remain sufficiently off-shell, and Eq. (j3J), with Eq. (JJJ) inserted, is free from 
the end-point singularity. 

The corresponding Sudakov factor can be derived in PQCD as a function of the transverse 
separation b and of the momentum fraction x carried by the spectator quark j^j, whose behavior 
is shown in Fig. EI The Sudakov factor suppresses the large b region, where the quark and the 
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Figure 3: QCD demands the presence of a Sudakov factor, which is an amplitude for a quark- 
antiquark color dipole not to emit real gluons in the final state. For a large transverse separation 
b, the quark and the antiquark do not shield each other's color charge, and intend to radiate. 
In this region Sudakov suppression is strong. 

antiquark are separated by a large transverse distance and the color shielding is not effective. 
It also suppresses the x ~ 1 region, where a quark carries all of the meson momentum, and 
intends to emit real gluons in hard scattering. The Sudakov factors from resummation [22] 
for the B and mesons are only associated with the light spectator quarks, since the double 
logarithms arise from the overlap of the soft and mass (collinear) divergences. These factors, 
being universal, are the same as in all our previous analyses. 

Similarly, the small x region corresponds to a configuration with a soft spectator, i.e., with 
a large color dipole in the longitudinal direction. The probability for this large color dipole not 
to radiate in hard scattering is also described by a Sudakov factor, which comes from threshold 
resummation for the hard kernels. For the derivation of this Sudakov factor, refer to [211] • For 
convenience, it has been parametrized as [21], 



with the constant c = 0.35. The above parametrization is motivated by the qualitative behavior 
of St' St(x) —>■ as x — > 0, 1 [211] ■ Since threshold resummation is associated with the hard 
kernels, the result could be process-dependent. It has been observed [23] that its effect is 
essential for factorizable decay topologies, and negligible for nonfactorizable decay topologies. 

2.3 Factorization Assumption 

We review the basics of FA for the B — > D^*'M decays. The relevant effective weak Hamiltonian 
is given by 





n eS = v cb v: d \cMOM + c 2 (//)o 2 (/i) 



(9) 



where the four-fermion operators are 





(10) 
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with the definition (q~xq 2 )v-A = 917^(1 — 75)92, V's the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements, and C\ and C 2 the Wilson coefficients. The B° — > D + n~ mode is referred to 
as the class- 1 (color-allowed) topology, in which the charged pion is emitted at the weak vertex. 
The B° — ► D°7T mode is referred to as the class-2 (color-suppressed) topology, in which the 
D° meson is directly produced. 

For the B° — > D + tt~ mode, 2 and Fierz transformed 0\ contribute. For the B° — > D°7r° 
mode, 0\ and Fierz transformed 2 contribute. Applying FA PQ or generalized FA (2112111]** 
to the hadronic matrix elements, we have 

(D+7r-\(cb)v-A(du) v ^ A \B°) « (D + \(cb) v _ A \B°)(n-\(du) v _ A \0) , 
(D\°\(db) v ^ A (cu) v „ A \B°) « (^K^v-aI^X^K^v-aIO) • (11) 

Substituting the definition of the B meson transition form factors, F BD and F Bw , and of the 
meson decay constants, f n and fo, the B° — ► £) + 7r~ ( class- 1) and 5° — > D°7r° (class-2) decay 
amplitudes are expressed as 

A(B°^D+n-) = t^V^Aml-mDUF^imDa^Dn), (12) 

V2^°-5V) = -&V*V: d {ml-ml)f D F B *{m 2 D )a 2 {Dit), (13) 

where the parameters a\ and 02 are defined by 

*x = CM + °&>, a 2 = C&) + °?jg>-, (14) 

N c being the number of colors. The B~ — > D°tt~ mode, involving both classes of amplitudes, 
is referred to as class-3. The isospin symmetry implies 

A(B° -> D+tT) = -> DV") + v/2 A(5° -> DV) . (15) 

It is straightforward to apply FA to other B — > D^M modes. 01 and a 2 depend on the color 
and Dirac structures of the operators, but otherwise are postulated to be universal [H I2H1 122 • 
They have the orders of magnitude ai(Dn) ~ 0(1) and a 2 (Di[) ~ 0(1/ N c ). The consistency 
of FA can be tested by comparing a\ and a 2 extracted from various decays. Within errors, the 
class- 1 decays B° — > D" + M~ with M = it, p, a 1; D s , and D* are described using a universal 
value \ai\ « 1.1 ±0.1, whereas the class-2 decays i? — > K^M with M = J/ip and ip(2S) suggest 
a nearly universal value \a 2 \ ~ 0.2-0.3 |28j . The wide range of \a 2 \ is due to the uncertainty 
in the B — > K^*' form factors. The class-3 decays B~ — > D^*'°M~ with M = ir and p, which 
are sensitive to the interference of the two decay topologies, can be explained by a real and 
positive ratio 02/01 ~ 0.2-0.3, which seemed to agree with the above determination of |oi| and 
1 0.2 1 - This is the reason FA was claimed to work well in explaining two-body charmed B meson 
decays, before the class-2 modes B° — ► D°M° with M = ir, 77, and u were measured. 

The recently observed B° — > D°M° branching ratios listed in Table ^ 123 120] revealed 
interesting QCD dynamics. The parameter \a 2 \ directly extracted from these modes falls into 

**The main difference between FA and generalized FA is that nonfactorizable contributions are included in 
the latter. 
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the range of |a2(-D7r)| ~ 0.35— 0.60 and |a2(-D*vr)| ~ 0.25—0.50 |3T| . To maintain the predictions 
for the class-3 decays, there must exist sizeable relative strong phases between class-1 and 
class-2 amplitudes [T5j . which are Arg(a 2 /ai) = 59° for the Dti modes and Arg(a>2/ai) = 63° 
for the D*ir modes [3*T] . These results can be regarded as a failure of FA: the parameters 
ct2 in different types of decays, such as B — > D^n and B — > K(*'J/ip, differ by almost a 
factor 2 in magnitude, implying strong nonuniversal nonfactorizable effects. It is then crucial 
to understand this nonuniversality and, especially, the mechanism responsible for the large 
relative phases in a systematic QCD framework. 



3 B — ► Dn IN PQCD 

In this section we take the B — > Dir decays as an example of the PQCD analysis. The intensive 
study of all other modes will be performed in the next section. The B — > Dir decay rates have 
the expressions, 

? l = ^Gl\V cb \*\V ud \^\A l \* . (16) 
1287T r 

The indices for the classes i = 1, 2, and 3, denote the modes B° — > D + n~, B° — > D°n°, and 
B~ — > D n~, respectively. The amplitudes Aj are written as 

M = f'Mxt + /iJ^cxc + -Mext + M cxc , (17) 

V2A 2 = -(f D U - IbUc + M int - M cxc ) , (18) 

^3 = fMxt + /o6nt + -Mext + -Mint , (!9) 

with Jb being the B meson decay constant. The functions £ ex t, £int, and ^ ex c denote the 
factorizable external iy-emission (color-allowed), internal J^-emission (color-suppressed), and 
^-exchange contributions, which come from Figs. 4(a) and 4(b), Figs. 5(a) and 5(b), and 
Figs. 6(a) and 6(b), respectively. The functions A4 ex t, M-mu and Ai cxc represent the nonfactor- 
izable external iy-emission, internal iy-emission, and VT-exchange contributions, which come 
from Figs. 4(c) and 4(d), Figs. 5(c) and 5(d), and Figs. 6(c) and 6(d), respectively. All the 
topologies, including the factorizable and nonfactorizable ones, have been taken into account. 
It is easy to find that Eqs. (fT7]) - (|T9"j) obey the isospin relation in Eq. (fTB]) . 




(a) (b) (c) (d) 



Figure 4: Color-allowed emission diagrams contributing to the B — > D^*'tt decays. 

In the PQCD framework based on factorization theorem, an amplitude is expressed as 
the convolution of hard b quark decay kernels with meson wave functions in both the longitu- 
dinal momentum fractions and the transverse momenta of partons. Our PQCD formulas are 
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(a) (b) (c) (d) 



Figure 5: Color-suppressed emission diagrams contributing to the B — > D^ir decays. 




(a) (b) (c) (d) 



Figure 6: Annihilation Diagrams contributing to the B — > D^tc decays. 

derived up to leading-order in a s , to leading power in m D /m B and in A/m D , and to lead- 
ing double-logarithm resummations. For the Wilson coefficients, we adopt the leading-order 
renormalizat ion-group evolution for consistency, although the next-to-leading-order ones are 
available in the literature [32] • For the similar reason, we employ the one- loop running cou- 
pling constant a s (fi) = 27r/[ / 9 1 hi(fi/ Aqq D )] with (5\ = (33 — 2nj)/3, rif being the number of 

(5) 

active quarks. The QCD scale is chosen as Aqq D = 193 MeV for the scale m b < jj, < mw, 

which is derived from Aqq D = 250 MeV for jj < m^. 

The leading-order and leading-power factorization formulas for the above decay amplitudes 
are collected in Appendix A. Here we mention only some key ingredients in the calculation. 
The formulas for the B —>■ decays turn out to be simpler than those for the B — ► titt ones. 
The simplicity is attributed to the power counting rules under the hierachy of the three scales in 
Eq. The hard kernels are evaluated up to the power corrections of order A/mo ~ m^/m^ 
(A/ tub is regarded as being of even higher power). Following these rules, the terms proportional 
to x\ ~ A/bib and to X2 ~ A/wtD are higher-power compared to the leading 0(1) terms and 
dropped. We have also dropped the terms of higher powers in r = m D /m B . Accordingly, 
the phase space factor 1 — r 2 , appearing in Eq. (|16|) originally, has been approximated by 1. 
This approximation is irrelevant for explaining the ratios of the B — > Dtt branching ratios, and 
causes an uncertainty in the absolute branching ratios, which is much smaller than those from 
the CKM matrix element \V C b\, and from the meson decay constants J'b and fp. 

Up to the power corrections of order A/ms and A/mj, we consider only a single B (D) 
meson wave function. The nonperturbative B meson and pion wave functions have been fixed 
in our previous works ^2JE3- The unknown D meson wave function was determined by fitting 
the PQCD predictions for the B — ► D transition form factors to the observed B — > Dlv decay 
spectrum [211] • The contributions from the two-parton twist-3 D meson wave functions, being 
higher-power, are negligible. Note that in the charmless decays the contributions from the 
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two-parton twist-3 light meson distribution amplitudes are not down by a power of I/tub- 
These distribution amplitudes, being constant at the momentum fraction x — > as required 
by equations of motion, lead to linear singularities in the collinear factorization formulas. The 
linear singularities modify the naive power counting, such that two-parton twist-3 contributions 
become leading-power [2*4^ 13*3*] . In the charmed decays the above equations of motion are 
modified [201 , and the two-parton twist-3 D meson wave functions vanish at the end point of 
x. Therefore, their contributions are indeed higher-power. 

Retaining the parton transverse momenta kx, the nonfactorizable topologies generate strong 
phases from non-pinched singularities of the hard kernels For example, the virtual quark 
propagator in Fig. 5(d) is written, in the principle- value prescription, as 



x 3 (x 2 - Xi)m 2 B - (k 2T - ki T + k 3T ) 2 + it 



P 



1 



x 3 (x 2 - xi)m B - (k 2T - kiT + k 3T ) 2 _ 
-in5(x 3 (x 2 - x 1 )m B - (k 2T - k 1T + k 3T ) 2 ) , 

(20) 

with £3 being the momentum fraction associated with the pion. The second term contributes 
to the strong phase, which is thus of short-distance origin and calculable. The first term in the 
above expression does not lead to an end-point singularity. Note that the strong phase from 
Eq. ({2D}) is obtained by keeping all terms in the denominator of a propagator without neglecting 
x\ and x 2 . 



4 NUMERICAL RESULTS 

The computation of the hard kernels in kj- factorization theorem for other charmed decay modes 
is similar and straightforward. The B —>■ D*tt decay amplitudes are the same as the B — > Dtx 
ones but with the substitution of the mass, the decay constant, and the distribution amplitude, 

m D ^m D *, f D ^fD*, <Pd(x 2 ) -> (j) D *{x 2 ) . (21) 

This simple substitution is expected at leading power under the hierachy in Eq. (0): the dif- 
ference between the two channels should occur only at 0(A/m_o). An explicit derivation shows 
that the difference occurs at the twist-3 level for the nonfactorizable emission diagrams in 
Fig. 4 and for the annihilation diagrams in Fig. 6. It implies that the universality (channel- 
independence) of at and a 2 assumed in FA breaks down at subleading power even within 
the B -> £)WM decays. 

Replacing the pion in Figs. 4-6 by the p (u) meson, we obtain the diagrams for the B — > 
D^p(u) decays. The factorization formulas for the B — > Dp(u>) decay amplitudes are also the 
same as those for the B — > Dtt ones but with the substitution, 

-> <P P) u , € ^ , <f>l -> # lW > ^0 -> m pH , (22) 

where <fi n (<$f) is the two-parton twist-2 (twist-3) pion distribution amplitude, (pp^ (0p^) the 
two-parton twist-2 (twist-3) p and u meson distribution amplitudes, respectively, mo the chiral 
enhancement scale, and m^y the p (u) meson mass. Hence, the similar isospin relation holds: 

A(B° -> D-p + ) = A(B + -> D°p + ) + V2A{B° -> D°p°) . (23) 
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The B — > D*p(uj) decays contain more amplitudes associated with the different polariza- 
tions. However, at leading power, the amplitudes associated with the transverse polarizations, 
suppressed by a power of m D */m B or of m p /m B , are negligible. That is, the factorization for- 
mulas for the B — > D*p(uj) modes are the same as the B — > Dp(u) ones with the substitution 
in Eq. (}2*Tj) . Note that the H^-exchange contribution changes sign in the B° -> DM°u decay 
amplitude: 

A (B° - D^°u) = ~(/ DW ^ + feUc + -Mint + -Mcxc) , (24) 

due to the different quark structures between the u meson (proportional to uu + dd) and the 
p° meson (proportional to uu — dd). 

In the numerical analysis we adopt the model for the B meson wave function, 



4> B (x, b) = N B x 2 (l — x) 2 exp 



H 


( xm B s 


2 , ,2 1.2" 




V cob s 


1 2 



(25) 



with the shape parameter uj b and the normalization constant iV^ being related to the decay 
constant f B through 

/**.<*.<>) = j^jj. (26) 

The Z)W meson distribution amplitude is given by 

= i/i)W*(l - x)[l + C DW (1 - 2x)] , (27) 

with the shape parameter C D (*). The pion and p(cu) meson distribution amplitudes have been 
derived in jSSllSZl, whose explicit expressions are given in Appendix A. The B meson wave 
function was then extracted from the light-cone-sum-rule (LCSR) results of the B — » 7r transi- 
tion form factor [23]. The range of C D {*) was determined from the measured B — > D^*'lv decay 
spectrum at large recoil employing the B meson wave function extracted above. We do not 
consider the variation of </>£>(*) with the impact parameter b, since the available data are not 
yet sufficiently precise to control this dependence. 
The input parameters are listed below: 

f B = 190 MeV , uj b = 0.4 GeV , 
f D = 240 MeV , C D = 0.8 ± 0.2 , 
f D , = 230 MeV , C D * = 0.7 ± 0.2 , 

U= 132 MeV, f p = f w = 200 MeV, f p = f w = 160 MeV , 

m B = 5.28 GeV , m b = 4.8 GeV , 

m D = 1.87 GeV , m D * = 2.01 GeV , m c = 1.3 GeV , 



m p = 0.77 GeV , m u = 0.78 GeV 



m t = 170 GeV , m w = 80.41 GeV , m = 1.4 GeV , 
r B ± = 1.674 x 10" 12 s , r B o = 1.542 x 10~ 12 s , 

G F = 1.16639 x 10" 5 GeV -2 , \V cb \ = 0.043 , \V ud \ = 0.974 , (2* 
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where m t , myy, Tb± and Tgo denote the top quark mass, the W boson mass, the B ± meson life- 
time, and the B° meson lifetime, respectively. The above meson wave functions and parameters 
correspond to the form factors at maximal recoil, 

F B7r ~0.3, ^ D ~0.57, ^- 0.52, (29) 

which are close to the results from QCD sum rules [HE1EH]- We stress that there is no arbitrary 
parameter in our calculation, though the value of each parameter is only known up to a range. 

The PQCD predictions for each term of the B — > Dn decay amplitudes are exhibited in 
Table [21 The theoretical uncertainty comes only from the variation of the shape parameter for 
the D meson distribution amplitude, 0.6 < Co < 1.0. It is expected that the color-allowed 
factorizable amplitude /^ext dominates, and that the color-suppressed factorizable contribution 
/spirit is smaller due to the Wilson coefficient Ci+C^/iVc ~ 0. The color-allowed nonfactorizable 
amplitude -M ex t is negligible: since the pion distribution amplitude is symmetric under the 
exchange of x% and 1 — 23, the contributions from the two diagrams Figs. 4(c) and 4(d) cancel 
each other in the dominant region with small X2- It is also down by the small Wilson coefficient 
Ci/N c . For the color-suppressed nonfactorizable contribution .Mint, the above cancellation does 
not exist in the dominant region with small 23, because the D meson distribution amplitude 
4>d{x2) is not symmetric under the exchange of X2 and 1 — X2. Furthermore, .Mint, proportional 
to C2/N c ~ 0.3, is not down by the Wilson coefficient. It is indeed comparable to the color- 
allowed factorizable amplitude f n £, ex t, and produces a large strong phase as explained in Eq. (j2T"j) . 
Both the factorizable and nonfactorizable annihilation contributions are small, consistent with 
our argument in Sec. II. 

The predicted branching ratios in Table El are in agreement with the averaged experimental 
data |2H1 HUB EDI- We extract the parameters a\ and 02 by equating Eqs. (""""2*j) and (JT3J) to 
Eqs. (II 7 j) and (["FT"") , respectively. That is, our a\ and 02 do not only contain the nonfactorizable 
amplitudes as in generalized FA, but the small annihilation amplitudes, which was first discussed 
in |"""T]. We obtain the ratio \a2/ax\ ~ 0.43 with 10% uncertainty and the phase of 02 relative 
to <2i about Arg(a2/ai) ~ —42°. If excluding the annihilation amplitudes /e^exc and -M ex c, we 
have 1 0,2/(11 1 ~ 0.46 and Arg(ci2/ 'ai) ~ —64°. Note that the experimental data do not fix the 
sign of the relative phases. The PQCD calculation indicates that Arg^/a^ should be located 
in the fourth quadrant. It is evident that the short-distance strong phase from the color- 
suppressed nonfactorizable amplitude is already sufficient to account for the isospin triangle 
formed by the B — > Dir modes. The conclusion that the data hint large final-state interaction 
was drawn from the analysis based on FA [T"*"l I3*T| 14*2*1 1^3*] . Hence, it is more reasonable to 
claim that the data just imply a large strong phase, but do not tell what mechanism generates 
this phase |44j . From the viewpoint of PQCD, this strong phase is of short distance, and 
produced from the non-pinched singularity of the hard kernel. Certainly, under the current 
experimental and theoretical uncertainties, there is still room for long-distance phases from 
final-state interaction. 

The PQCD predictions for the B — > D*n decay amplitudes and branching ratios in Tabled 
are also consistent with the data j""f"""J. Since m^* and are only slightly different from 
m D and (f> D , respectively, the results are close to those in Table El The B — > D*n branching 
ratios are smaller than the B — > Dir ones because of the form factors ^ff* < £+ D as shown in 
Eq. (J29J) . Similarly, the ratio (02/01 1 and the relative phase Ar g (c^/ai) are also close to those 
associated with the B — > Dtt decays. We obtain the ratio |a2/oi| ~ 0.47 with 10% uncertainty 
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and the relative phase about Arg(a 2 /ai) ~ —41°. Excluding the annihilation amplitudes, we 
have |a 2 /ai| ~ 0.5 and Arg(a 2 /ai) ~ —63°. 

The PQCD predictions for the B — > p{ui) branching ratios are listed in Table El which 
match the data jl^j. The B° — > D + p~ and B~ —>■ D°p~ branching ratios are about twice of the 
B° — > D + n~ and B~ — > D°ti~ ones because of the larger p meson decay constant, (fp/fn) 2 ~ 2. 
The relatively smaller B° — > D°p° branching ratio is attributed to the cancellation of the above 
enhancing effect between the color-suppressed and jy-exchange contributions, consistent with 
the observation made in an analysis based on the topological-amplitude parametrization |1d] . 
The B° — > D u branching ratio is larger than the B° — > D°p° one due to the constructive 
inteference between the color-suppressed contribution and the annihilation contribution as in- 
dicated in Eq. (}24j) . To obtain the B° — > D*p helicity amplitudes and their relative phases 
|Hj, the power-suppressed contributions from the transverse polarizations must be included. 
For consistency, the contribution from the longitudinal polarization should be calculated up to 
the same power. We shall study this subject in a forthcoming paper. The predictions for the 
B — > D*°p°, D*°uj decays can be compared with future measurement. The latter branching 
ratio is larger than the former one because of the same reason as for the B — > D°p°, D°u 
decays. 



5 COMPARISON WITH OTHER APPROACHES 

In this section we make a brief comparison of the PQCD formalism with other QCD approaches 
to exclusive B meson decays, emphasizing the differences. For more details, refer to 48]. As 
mentioned before, there are two kinds of factorization theorem for QCD processes [T^j: collinear 
factorization, on which the soft-collinear effective theory (SCET) |4UJ ED] , LCSR j^H E2] , and 
QCDF j^j are based, and kx factorization, on which the PQCD approach is based. Calculating 
the B — ► 7i form factor F Bn in collinear factorization up to leading power in 1/wig and leading- 
order in a s , an end-point singularity occurs. Hence, we define three types of contributions: a 
genuine soft contribution / s , a contribution / EP containing the end-point singularity, and a 
finite contribution / F : 

FS7r = /S + / EP + / F (3Q) 

The second term can not cover the complete soft contribution, because it is from a leading 
formalism. Note that the end-point singularity exists even in the heavy quark limit. Hence, 
B meson decays differ from other exclusive processes, which become calculable in collinear 
factorization at sufficiently large momentum transfer. 

There are two options to handle the above end-point singularity first, an end-point 
singularity in collinear factorization implies that exclusive B meson decays are dominated by 
soft dynamics. Therefore, a heavy-to-light form factor is not calculable, and f EP should be 
treated as a soft object, like / s . In SCET and QCDF, F B?r is then written, up to 0(a s ), as 

[Ml ESI 

F m = f F + f , (31) 

with 

/NF = /S + / EP ; 

f = B ®T'®0, . (32) 
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The soft form factor / NF , obeying the large-energy symmetry relations [56j, can be estimated 
in terms of a triangle diagram without a hard gluon exchange in LCSR j^BJEj. However, since 
the pion vertex has been replaced by the pion distribution amplitudes under twist expansion, 
what is calculated in LCSR is not the full soft contribution. The term / F has been expressed 
as a convolution of the hard-scattering kernel T' with the light-cone distribution amplitudes 
of the B meson and of the pion in the momentum fractions, implying that it is calculable in 
collinear factorization. 

The second option is that an end-point singularity indicates the breakdown of the collinear 
factorization. Hence, the kx factorization is the more appropriate framework, in which the 
parton transverse momenta kx are retained in the hard kernel, and / EP does not develop an 
end-point singularity. Both / EP and / F are then calculable, and expressed, in the PQCD 
approach, as 

F fl ' r = / EP + / F = ^®T®^, (33) 

where the symbol <8> represents the convolution not only in the momentum fractions, but in the 
transverse separations. The hard kernel T' in Eq. ()32)1 is derived from the complete hard kernel 
T by dropping the terms which lead to the end-point singularity in the collinear factorization. 
Certainly, the subtraction of these terms depends on a regularization scheme [S3]. The strong 
Sudakov suppression in the soft parton region implies that the genuine soft contribution / s 
is not important (HI 121] • Equation is then claimed to be a consequence of the hard- 
dominance picture, because a big portion of F Bn is calculable. The agreement between the 
sum-rule and PQCD predictions for many B meson transition form factors justifies that / s 
is indeed negligible. Since / EP remains in Eq. (|3*3*|) . the form factor symmetry relations at 
large recoil are still respected in the PQCD framework [21], which are then modified by the 
subleading term / F . 

Therefore, the soft-dominance (hard-dominance) picture postulated in LCSR (PQCD) makes 
sense in the collinear {kx) factorization [IB]- The two pictures arise from the different theoret- 
ical frameworks, and there is no conflict at all. In other words, the soft contribution refers to 
/ NF in SCET, LCSR, and QCDF, which is large, but to / s in PQCD, which is small. LCSR 
can be regarded as a method to evaluate / EP or / NF in the collinear factorization (at least the 
light-cone distribution amplitudes have been employed on the pion side), while the hp factor- 
ization is adopted in PQCD for the evaluation of / EP . We emphasize that there is no preference 
between the two options for semileptonic B meson decays, both of which give similar results as 
stated above. However, in their extension to two-body nonleptonic B meson decays, predictions 
could be very different. For example, the main source of strong phases in the B — > tttt decays 
is the correction to the weak vertex in QCDF, but the annihilation diagram in PQCD. This is 
the reason QCDF (PQCD) predicts a smaller and positive (larger and negative) CP asymmetry 
CVtt fBl EH] • It is then possible to discriminate experimentally which theoretical framework 
works better. 

Next we compare our formalism for two-body charmed nonleptonic B meson decays based 
on the kp factorization with SCET and QCDF based on the collinear factorization. Currently, 
LCSR has not yet been applied to the nonfactorizable contribution discussed here, but only 
to that from three-parton distribution amplitudes [21], since the former, involving two loops, 
is more complicated to analyze. Similarly, the neglect of kx results in end-point singularities 
in the factorizable contributions £ ex t and £i nt , which need to be parametrized in terms of the 
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B — > D and B — > tt transition form factors, respectively. It also causes an end-point singularity 
in the color-suppressed nonfactorizable amplitude .Mint, if the c quark is treated as being 
massive. This is why the color-suppressed modes, i.e., the magnitude and the phase of a 2 , can 
not be predicted in QCDF, and the proof of QCDF in the SCET formalism 60J considered 
only the color-allowed mode B° — > D + n~ . The color-allowed nonfactorizable amplitude A4 ex t 
is calculable in QCDF, because the end-point singularities cancel between the pair of diagrams, 
Figs. 4(c) and 4(d). We mention a recent work on SCET [UT], in which the color-suppressed 
nonfactorizable amplitude has been parametrized as an expression similar to Eq. (J31|) . 

If the c quark is treated as being massless, the end-point singularities in the pair of color- 
suppressed nonfactorizable diagrams, Figs. 5(c) and 5(d), will cancel each other as in the 
charmless case |38{ 15 8|. This can be understood by examining the behavior of the integrand of 
.Mint i n Eq. (|45|) in the dominant region with small 23, noticing that the D meson distribution 
amplitude 0d(^2) would be symmetric under the exchange of Xi and 1 — x 2 in the m c — > 
limit. However, the nonfactorizable contribution will become negligible in this limit, such that 
the amplitude .Mint, though calculable in QCDF, is not large enough to explain the B — > Dtx 
data. It is then obvious that the PQCD approach has made a great contribution here: the 
nonfactorizable corrections to the naive factorizations of both the color-allowed and color- 
suppressed modes can be predicted, and the latter is found to be very important. 

6 CONCLUSION 

In this paper we have analyzed the two-body charmed nonleptonic decays B — > D^*'M with 
M = it, p, and u in the PQCD approach. This framework is based on kx factorization theorem, 
which is free of end-point singularities and gauge-invariant, factorization theorem is more 
appropriate, when the end-point region of a momentum fraction is important, and collnear 
factorization theorem breaks down. By including the transverse degrees of freedom of partons 
in the evaluation of a hard kernel, and the Sudakov factors from and threshold resummations, 
the virtual particles remain sufficiently off-shell, and the end-point singularities do not exist. 
We have explained that there is no conflict between LCSR with the soft- dominance picture 
and PQCD with the hard- dominance picture, since the soft contributions refer to the different 
quantities in the two theoretical frameworks. 

The derivation of the factorization formulas for the B — ► D^M decay amplitudes follow 
the power counting rules constructed in our previous work on the B — > transition form 
factors. Under the hierachy mg ^> m D (*) ^> A, the B and meson wave functions exhibit 
a peak at the momentum fractions around A/m# and A/m D (*), respectively. Up to leading 
power in m D (*)/rriB and in A/m D (*), only a single B meson wave function and a single D'*' 
meson wave function are involved. The factorization formulas then become simpler than those 
for the charmless decays. Moreover, the factorization formulas for all the B — > D^*'M modes 
are identical, except the appropriate substitution of the masses, the decay constants, and the 
meson distribution amplitudes. We emphasize that there is no arbitrary parameter in our 
analysis (there are in QCDF), though all universal inputs are not yet known precisely. The 
meson wave functions have been determined either from the semileptonic data or from LCSR. 

Being free from the end-point singularities, all topologies of decay amplitudes are calculable 
in PQCD, including the color-suppressed nonfactorizable one. This amplitude can not be 
computed in QCDF based on the collinear factorization theorem due to the existence of the 
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end-point singularities for a massive c quark. We have observed in PQCD that this amplitude, 
not suppressed by the Wilson coefficient (proportional to C 2 /iV c ), is comparable to the dominant 
color-allowed factorizable amplitude. It generates a large strong phase from the non-pinched 
singularity of the hard kernel, which is crucial for explaining the observed B — > D^M branching 
ratios. The other topologies are less important: the color-allowed nonfactorizable contribution 
is negligible because of the pair cancellation and the small Wilson coefficient C\/N c . The color- 
suppressed factorizable amplitude with the small Wilson coefficient a 2 = C\ + C<ijN c is also 
negligible. The annihilation amplitudes are small, since they come from the tree operators. 

All our predictions are consistent with the existing measurements. For those without data, 
such as the B — > D*°p°, D*°u> modes, our predictions can be confronted with future measure- 
ment. As stated before, we have predicted the large strong phases from the scalar-penguin 
annihilation amplitudes, which are required by the large CP asymmetries observed in two-body 
charmless decays. The success in predicting the storng phases from the color-suppressed non- 
factorizable amplitudes for the two-body charmed decays further supports the kx factorization 
theorem. The conclusion drawn in this work is that the short-distance strong phase is already 
sufficient to account for the B — > D^*'M data. Certainly, there is still room for long-distance 
strong phases from final-state interaction. For the application of the PQCD approach to other 
charmed decays, such as B — > D^K and B — ► £H*)/ , refer to (52] and (53], respectively. 
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A FACTORIZATION FORMULAS FOR B —> Dn 

In this Appendix we present the factorization formulas for the B — > Dn decay amplitudes. We 
choose the B meson, D meson, and pion momenta in the light-cone coordinate as, 

P 1 = ^(1,1,0 T ), P 2 = ^(l,r 2 ,0 T ), P 3 = ^f(0,l-r 2 ,0 T ), (34) 

respectively, with r = m^/m^ being defined before. The fractional momenta of the light 
valence quarks in the B meson, D meson and the pion are 



fa = 0*^(1,0, Or) +k 1T for £ int , M- mt , 



fa = xi— ^(0, 1, Or) + kir for others, 

m B 

k 2 = ^ 2 -^|(l,0, Or) +k 2T , 
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k 3 = x 3 ^(0,l-r 2 ,0 T ) + k 3 T 



(35) 



respectively. Which longitudinal component of k\, kf or , is relevant depends on the final- 
state meson the hard gluon attaches. That is, it is selcted by the inner product k\ ■ k 3 or 
ki ■ k 2 . 

The factorizable amplitudes £ ext , £mt and £ exc are written as 



rl rl/A 

= 167iCFy/rrn 2 B / dxidx 2 / bidbib 2 db 2 <p B (xi,bi)(f)D(x 2 ) 
Jo Jo 

x E e (t<p)h(x u x 2 , bx, b 2 )S t {x 2 ) + rE e {tf)h{x 2 , x u b 2 , b^S^xx) 

rl rl/A 

£ int = IQ-kC F y/rm 2 B / dxidx 3 / bidbib 3 db 3 4> B (xi, bi) 

Jo Jo 

x{[(l + x 3 )Mx 3 ) + r„(l - 2x3)^(2:3) + KM)] 
xE i {tf ) )h(x 1 ,x 3 {l - r 2 ), 61, 6 3 )5 t (x 3 ) 
+2r o 0£(x 3 )^(tf } )% 3 , - r 2 ), 63, , 

rl rl/A 

^exc = l6nC F y/rm B / dx 2 dx 3 / b 2 db 2 b 3 db 3 (p D (x 2 ) 

Jo Jo 

x -x 3 7r (x 3 )£; a (4 1) )/i o (x2, x 3 (l - r 2 ), 6 2 , ^^(a^) 
+x 2 7r (a; 3 ) J E a (^ 1) )/i a (a; 3 , x 2 (l - r 2 ), 6 3 , 6 a )S , t (a;a) 

with the mass ratio r = m /m B , the evolution factors, 

£ e (t) = a s (t) ai (t)ex P [-S B (t)- S D (t)] , 
Ei(t) = a s (t)a 2 (t) exp[-S B (t) - S n (t)] , 
£ (t) = a s (t)a 2 (t)exp[-S D (t) - S„(t)] , 



(36) 



(37) 



(3* 



and the Wilson coefficients, 



ai = C 2 



«2 = Ci + 



C2 



(39) 



(40) 



Note that C\ = and C 2 = 1 at tree level in our convention. The explicit expressions of the 
Sudakov factors exp[— S B (t)}, exp[— S^t)] and exp[— S n (t)] from fcy resummation are referred 

to [2011211. 

The functions h's, obtained from Figs. 4(a) and 4(b), Figs. 5(a) and 5(b), and Figs. 6(a) 
and 6(b), are given by 



h(x 1 ,x 2 ,b 1 ,b 2 ) 



h a (x 2 ,x 3 ,b 2l b 3 ) 



K Q {y/x 1 x 2 m B bx) 

x [6{bi - b 2 )K (^fx~ 2 m B bx) I {^fx~ 2 m B b 2 ) 
+9(b 2 - bx)Ko (y/x^m B b 2 ) I {y/xim B bi)] 



(41) 



i 7 ^) (^x 2 x 3 m B b 2 ) 



9{b 2 - b 3 )H^ ] (y/x^m B h) J (yfx 3 ~m B b 3 
+6(b 3 - b 2 )H^ ] (yfx3~m B b 3 ) J {y/xim B bz) 



(42) 
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The hard scales t are chosen as 

t)}> = max( v /i 2 ~m jB , l/&i, l/6 2 ) , t^ 2) = max(^/x^m B , 1/bx, l/6 2 ) , 

tf } = max(^x 3 (l - r 2 )m B , l/&i, 1/63) , 4 2) = max(y / xi(l - r 2 )m B , I/61, 1/&3) , 

= max( ^3(1 -r 2 )m B , l/6 2 , I/63) , if = max(^r 2 (l - r 2 )m B , l/b 2 , 1/63) . (43) 

For the nonfactorizable amplitudes, the factorization formulas involve the kinematic vari- 
ables of all the three mesons. Their expressions are 



M 



cxt 



A* 



hit 



rl rl/A 

32nV2NC F y/rm B / [dx] / bidbib 3 db 3 (f) B (xi ) bi)<i) D (x 2 )(j) 1x (x 3 ) 
Jo Jo 

x [ar 3 J5 6 (4 1) )/i6 1) (a:i,6i) - (1 - ^3 + x 2 )E b {tf ) )hf ) {x h , 

rl rl/A 

327rv 2NCF\frm B \ [dx] / 61 d6i 6 2 d6 2 B (xi, 61 )</>£> (:r 2 ) 
Jo Jo 

x (-x 2 - x 3 )0 w (x 3 ) J E; d (4 1) )4 1) (x i ,6 i ) + (1 - x 2 )(p^x 3 )E d {tf)hf{x l ,b i 

rl rl/A 

32ixv2NCF\/rm B / [dx] / bidbib 2 db 2 (f) B (x 1 ,bi)(f)£)(x 2 ) 
Jo Jo 

x 3 (j)Tr(3c 3 )Ef(tf')hy'(xi,bi) - x 2 <fin(x 3 )Ef(tf>)hf\xi,bi 



(44) 



,(45) 



(46) 



from Figs. 4(c) and 4(d), Figs. 5(c) and 5(d), and Figs. 6(c) and 6(d), respectively, with the 
definition [dx] = dx\dx 2 dx 3 . The evolution factors are given by 

E b (t) = ^ s (t)^-exp[-S(t)\ b2=bl ], 

E d (t) = a s (t)^exp[-S(t)\ b3=bl ], 

E f (t) = a s (t)^exp[-S(t)\ h=b2 ]. 

with the Sudakov exponent S = S B + S B + S n . 

The functions j — 1 and 2, appearing in Eqs. (|44j) - (J46j) . are written as 



(47) 



h 



h 



00 



[0(6i - 63)^0 (Sm fl 6i) Jo (Bm B b 3 ) 
+6(b 3 - bi)K (Bm B b 3 ) I Q (Brush)] 
( K {B jmB b 3 ) for B) > \ 

X ^ fH^(^\Bj\m B b 3 ) for 5 2 < J 

[eih-b^KoiDmMlaiDmM 
+6(b 2 - h)K (Dm B b 2 ) J pm B &0] 
/ K (D jmB b 2 ) for £>? > \ 
fF«( A /|D?|m B 6 2 ) for D 2 < J 



(46 



x 

7T 



(49) 



2 "U \y |- 3 i"»B"zy "J 

0(6! - 6 2 )tf ;i) (Fm B bi) J (Fm B b 2 ) 
-0(&2 - 6i)^ 1} (Fm B b 2 ) J (Fm B bi) 



ITT TT 



(i) 



for F 2 > 



|F. 2 |m B 6i) for F 2 < 



(50) 
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with the variables 



n 2 

B 


— X\X 2 , 




B{ 


= X\X 2 - X 2 X 3 {1 


-r z ) , 


Bl 


= X\X 2 - X 2 (l - 


x 3 )(l-r 2 ) , 


D 2 


= Xix 3 (l - r 2 ) , 




Dl 


= F 2 = (xi - x 2 


)x 3 (l-r 2 ) , 


Dl 


= (xi + x 2 )r 2 - 


(1 - x x - a; 2 )ar 3 (l - r 2 ) , 


F 2 


= x 2 x 3 (l - r 2 ) , 




F 2 


= X\ + x 2 + (1 - 


X\ - x 2 )x 3 (l - r 2 ) . 



(51) 

There is an ambiguity in defining a light-cone B meson wave function for the nonfactorizable 
amplitude -M e xc, since both the components kf and contribute through the inner products 
k\ ■ k 2 and k\ ■ fc 3 in the denominators of the virtual particle propagators. However, a careful 
examination of the factorization formula shows that the dominant region is the one with k 2 ~ 
0(A) and /c 3 ~ O(m^) at leading twist. Hence, we drop the term k\ ■ k 2 . The scales are 
chosen as 



4 J) = max(Bm B , ^J\B 2 \m B , 1/6 1; I/63) 
4 J) = max(Dm B , J\D 2 \m B , 1/& X , l/b 2 ) 



max{Fm B , J\F 2 \m B , l/b u l/b 2 ) . (52) 



We explain that the factorization formulas presented above are indeed of leading power under 
the power counting rules in [20]. The factorizable amplitudes are as shown in fK)\ I33j . For the 
nonfactorizable amplitudes, the terms proportional to x 3 and to 1 — x 3 in A4 ext cancel each 
other roughly. This cancellation can be understood by means of the corresponding expression 
in collinear factorization theorem: the first and second terms in Ai ex t are proportional to 

£3 l-x 3 + x 2 

(53) 



x x x\x 3 ' x x x\{\ - x 3 ) ' 

For simplicity, x\ has been suppressed, when it appears in the sum together with x 2 or x 3 . It 
is found that the first ratio cancels the 1 — x 3 term in the second ratio. That is, the x 2 term 
is in fact leading and not negligible. For a similar reason, the — x 2 term in cancels the 

1 — x 2 term. Hence, the — x 3 term is leading. If one drops — x 2 in .Mint, the above cancellation 
disappears, and a fake leading term will be introduced. 

The pion and p meson distribution amplitudes have been derived in EZ] : 

^(x) = -^h=x(l-x) [l + 0.44C 2 3/2 (2x-l) + 0.25C 4 3/2 (2x-l)l , (54) 



V^Nc 

fi(x) = ^=[l + 0A3C 1 2 /2 (2x-l) + 0.09Cl /2 (2x-l)\ , (55) 
4(ar) = ^-^(l-2a;)[l + 0.55(10x 2 -10x + l)l , (56) 
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(f)p(x) 

<t> v p {x) 



^=x(l-x) [l + 0.18C 2 3/2 (2x-l) 



J o 



{3(2x - l) 2 + 0.3(2:r - l) 2 [5(2:r - l) 2 - 3] 



2v/2lVe 

+0.21[3 - 30(2:r - l) 2 + 35(2a; - l) 4 ]} , 
3f T 

p (l-2x) 1 + 0.76(10x 2 - 10x + 1) 



3/J 



x(l-x) [l + 0.2C 2 3/2 (2:r- 1) 



■[1 + (2x - l) 2 ] + 0.24[3(2:r - l) 2 - 1] 



4v/2iV; 

with the Gegenbauer polynomials, 



2y/W c U l 

+0.12[3 - 30(2x - l) 2 + 35(2x - 1) 4 ]| , 

3 f r 

J 9 (l-2x) 1 + 0.93(10x 2 - 10x + 1) 



C 2 1/2 W = |(3t 2 - 1) , Cl /2 (t) = l(35t 4 - 30t 2 + 3) 



cTH) = l^t 1 - 1) , Ct'\t) 



,3/2, 
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(21t 4 - 14t 2 + 1) 



(57) 

(58) 
(59) 

(60) 

(61) 
(62) 



(63) 



We shall assume that the uj meson wave functions are identical to the p meson ones in this 
work. 
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Decay mode 



Belle 29' 



CLEO 30^ 



TTW 
-0.32 

2.20^;1 ± 0.79 



B° 
B° 
B° 
B° 
B° 
B° 



D°7l° 

D*°n° 

D°T] 
D*°T] 
D°LU 
D*°LU 



3.1 ±0.4 ±0.5 

9 7+0.8+0.5 
z - ' -0.7-0.6 

1.4^ ±0.3 



2.0 



h0.9 
-0.8 



±0.4 



1.8 ±0.5 

1+1.3 



-0.4 
-0.3 



3.1Ii;i ±0.8 



2.74] 



±0.55 



Table 1: Data (in units of 10~ 4 ) of the B° -> L>(*)°M° (X = 7i,r],u) branching ratios. 



Amplitudes 


C D = 0.6 


C D = 0.8 


C D = 1.0 




6.90 


7.46 


8.01 




-1.44 


-1.44 


-1.44 




-0.01 — 0.03a 


-0.02 - 0.03i 


-0.02 - 0.03i 


•Mext 


-0.24±0.57i 


-0.25±0.60i 


-0.27±0.65i 


Mint 


3.34- 3.02i 


3.22 - 3.07z 


3.10 - 3.12i 


M cxc 


-0.26 - 0.89z 


-0.31-0.95i 


-0.37- 1.02i 



Table 2: Predicted B — > Dn decay amplitudes in units of 10 2 GeV. 
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Quantities 


C D = 0.6 


C D = 0.8 


C D = 1.0 


Data 


A 1 


6.39 - 0.35i 


6.88 - 0.38i 


7.35 - OAOi 




A 2 


—1.53 + 1.48i 


-1.49 + 1.48i 


-1.45 + 1.45i 




A 3 


8.56 - 2A5i 


8.99 - 2.47i 


9.40 - 2AQi 




B(B° -> D + n-) 


2.37 


2.74 


3.13 


3.0 ±0.4 


B(B° -> D°7T ) 


0.26 


0.25 


0.24 


0.29 ±0.05 




4.96 


5.43 


5.91 


5.3 ±0.5 


a 2 /ai (w/o anni.) 


0.47(0.51) 


0.43(0.46) 


0.39(0.42) 




Arg(a,2/ai) (w/o anni.) 


-42.5°(-61.5° ) 


-41.6° (-63.5°) 


-41.9°(-65.3°) 





Table 3: Predicted B — > Dft decay amplitudes in units of 10 2 GeV, branching ratios in units 
of 10~ 3 , \ci2/ai\, and relative angle Arg{a 2 j ' a.\) in units of degree. 



Quantities 


C D * = 0.5 


C D * = 0.7 


C D * = 0.9 


Data 


A 1 


6.32 - 0A2i 


6.81 - 0A5i 


7.30 - 0A9i 




A 2 


—1.65 + 1.61i 


—1.62 + 1.59z 


—1.59 + 1.57i 




A, 


8.65 - 2.69i 


9.10-2.70i 


9.55 - 2.70i 




B(B° -> D* + ir~) 


2.16 


2.51 


2.88 


2.76 ±0.21 


B (B° -> D*°n°) 


0.29 


0.28 


0.27 


0.25 ±0.07 


B(B~ -> D*°7T-) 


4.79 


5.26 


5.75 


4.60 ± 0.40 


02/01] (w/o anni.) 


0.52 (0.55) 


0.47 (0.50) 


0.43 (0.47) 




Arg(a 2 /oi) (w/o anni.) 


-40.5° (-61.4°) 


-40.7° (-63.1°) 


-40.8° (-64.8°) 





Table 4: Predicted B — > D*n decay amplitudes in units of 10 2 GeV, branching ratios in units 
of 10~ 3 , (02/01 1, and relative angle Arg(a 2 /ai) in units of degree. 



Branching ratios 


= 0.6 


C D = 0.8 


C D = 1.0 


Data 


B{B° -> D+p-) 


4.10 


4.72 


5.38 


8.0 ± 1.4 


B (B° _> D°p°) 


0.17 


0.17 


0.17 


0.29 ±0.11 


B(B~ -> D°p-) 


7.26 


8.15 


9.09 


13.4 ± 1.8 


B(B° -> DV>) 


0.50 


0.54 


0.56 


0.18 ±0.06 


Branching ratios 


C d* = 0.5 


C D * =0.7 


C D * = 0.9 


Data 


B(B° -> £>* + p~) 


5.32 


6.16 


7.08 


7.3 ± 1.5 




0.18 


0.18 


0.19 


< 0.51 


-> L>*V) 


9.14 


10.32 


11.60 


15.5 ±3.1 


5(5° -> L>*°w) 


0.49 


0.53 


0.58 


< 0.74 



Table 5: Predicted B — > p(uo) branching ratios in units of 10 3 . 
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